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Abstract. For infinite (bulk) quantum fluids of particles interacting via pairwise 
sufficiently smooth interactions, the Wigner-Kirkwood formalism provides a 
semiclassical expansion of the Boltzmann density in configuration space in even powers 
of the thermal de Broglie wavelength A. This result permits one to generate an 
analogous A-expansion for the bulk free energy and many-body densities. The present 
paper brings a technically nontrivial generalization of the Wigner-Kirkwood technique 
to semi-infinite quantum fluids, constrained by a plane hard wall impenetrable to 
particles. In contrast to the bulk case, the resulting Boltzmann density involves also 
position-dependent terms of type exp(— 2x 2 /X 2 ) (x denotes the distance from the wall 
boundary) which are non-analytic functions of the de Broglie wavelength A. Under 
some condition, the analyticity in A is restored by integrating the Boltzmann density 
over configuration space; however, in contrast to the bulk free energy, the semiclassical 
expansion of the surface part of the free energy (surface tension) contains odd powers 
of A, too. Explicit expressions for the leading quantum corrections in the presence of 
the boundary are given for the one-body and two-body densities. As model systems for 
explicit calculations, we use Coulomb fluids, in particular the one-component plasma 
defined in the i/-dimensional (integer v > 2) space. 
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1. Introduction 

In this paper, we consider a quantum system of N identical particles j = 1,2, ... ,N 
of mass m which position vectors r 1; r 2 , . . . , r^r are confined to a region A of the v- 
dimensional space R v . For the sake of brevity, we denote the i/iV-dimensional position 
vector in configuration space by r = (ri, r 2 , . . . , rjv) and the corresponding gradient by 
V = (Vi, V2, • • • , Vjv). The Hamiltonian of the particles (in the absence of a magnetic 
field) is 

#=^(-i^) 2 + nr), (1.1) 

where h is Planck's constant and the scalar potential V(r) includes pairwise particle 
interactions J2j<k v {\ r j ~ r k\) plus perhaps one-particle interactions J2j u { r j) with some 
external potential u(r). The equilibrium statistical mechanics of the particle system is 
studied in the canonical ensemble at the temperature T (or, alternatively, the inverse 
temperature j3 = with k-g being Boltzmann's constant). In the classical regime, 

the Hamiltonian (11. ip is replaced by H = p 2 / (2m)+V(r), where p = (pi, p 2 , . . . , Pat) is 
the vector in the classical momentum space. The integration of the Boltzmann density 
exp(—j3H) over the momentum variables results in a trivial constant, so that the classical 
Boltzmann density in configuration r-space oc exp[— j3V{f)]. Quantum effects are related 
to the thermal de Broglie wavelength A = h(/3 /m) 1 / 2 . A nearly classical regime of present 
interest is defined as being such that the dimensionless parameter X/l (I being a typical 
microscopic length of the classical particle system, proportional to the mean interparticle 
distance) is sufficiently small. 

For a nearly classical bulk system in an infinite //-dimensional space A = R u , the 
Boltzmann density in configuration space r can be expanded in powers of A 2 (odd powers 
of A vanish) within the well-known Wigner-Kirkwood expansion [TJ [2] . The coefficients 
of the A 2 -expansion of the partition function, which is given by the integration of the 
Boltzmann density over configuration space r, are expressible in terms of the classical 
averages of gradient operations applied to V(r) and of their products. 

The applicability of the Wigner-Kirkwood expansion scheme depends on the form 
of the interaction potential v(r). For certain potentials, the expansion in A 2 is analytic in 
the sense of holomorphy. There are potentials for which the Wigner-Kirkwood expansion 
is a limited Taylor expansion in A 2 : divergence of an expansion coefficient indicates a 
loss of analyticity of the corresponding and higher-order quantum corrections. Like for 
instance, for the 3D exponential, screened Coulomb and square barrier potentials, the 
leading expansion term was found to have the non-analytic form -v/A 2 [3] . In the case of 
an inverse-power-law repulsive potential v(r) = t> (a/r) ra with the exponent n from the 
range 1 < n < 00, the Wigner-Kirkwood expansion turns out to be analytic in A 2 in the 
holomorphic sense [3]. In the hard-core limit n — > 00, the Wigner-Kirkwood expansion 
fails and one has once more the nonanalyticity of type x/A 2 , as was shown in numerous 
analytic studies (see, e.g., [H El US E] ) ■ 
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Much effort has been devoted to the calculation of quantum corrections for 3D 
Coulomb fluids of ±e charges, in the high-temperature regime described adequately by 
the classical Debye-Hiickel theory. Quantum corrections of the bulk equation of state 
were derived and analyzed in Ref. pj. As concerns the two-body density, its leading 
quantum correction exhibits at large distances an inverse-power-law decay in contrast 
with the classical Debye-Hiickel exponential screening (see the review [9]). 

To our knowledge, there does not exist in the literature a boundary version of 
the Wigner-Kirkwood expansion applicable to fluids in domains A surrounded by hard 
walls impenetrable to particles (or, equivalently, there is an infinite external potential 
outside the domain A). Surface properties of quantum particle systems have been 
investigated only in the presence of some smooth confining external-potential barriers 
(see, e.g., [101 El)- T ne so ^ e exceptions are represented by semiclassical studies of 3D 
quantum Coulomb fluids in the high-temperature Debye-Hiickel limit for which one has 
the explicit results for the density profile near a plain hard wall [12] and for the large- 
distance tail of charge correlations along a plain hard wall or a conducting wall |13j . In 
these explicit results, quantum corrections linear in A together with non-analytic terms 
of type exp(— 2x 2 /A 2 ) (x denotes the distance from the boundary) appear; the latter are 
linked to the propagator of a free particle [H] . 

The present paper brings a technically nontrivial generalization of the Wigner- 
Kirkwood technique to semi-infinite quantum fluids, constrained by a plane hard wall. 
In contrast to the bulk case, the resulting Boltzmann density also involves position- 
dependent terms which are non-analytic functions of the de Broglie thermal wavelength 
A. Under some condition about the classical density profile, the analyticity in A is 
restored by integrating the Boltzmann density over configuration space; however, in 
contrast to the bulk free energy, the semiclassical expansion of the surface part of the 
free energy (surface tension) contains odd powers of A, too. Explicit expressions for the 
leading quantum corrections in the presence of the boundary are given for one-body and 
two-body densities. 

As a model system for explicit calculations, we shall use a special kind of Coulomb 
fluid, namely the one-component plasma (jellium) made of mobile pointlike charges e 
neutralized by a uniform oppositely charged fixed background, in the //-dimensional 
space (because of some physical nontrivialities in ID, we shall restrict ourselves to 
dimensions v > 2). In the infinite space A = R u , or when the hard walls surrounding the 
confining domain A do not induce image charges (plain hard walls), the translationally 
invariant Coulomb interaction potential of two e-charges v(r) with r = |r 3 - — r^| is defined 
as the solution of the Poisson equation V 2 t> = — s u e 2 5(r), where s u = 2tt ij ' 2 /T(u/2) (T 
stands for the Gamma function) is the surface area of the //-dimensional unit sphere. 
The resulting Coulomb interaction potential is logarithmic in v = 2 dimensions and of 
type e 2 r 2 ~7(z/-2) in v > 3 dimensions. The present definition of the Coulomb potential 
maintains many generic properties of "real" 3D Coulomb systems, like screening and the 
corresponding sum rules (see the review [15J). The total potential V(r) then satisfies 
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for each of the particle coordinates the differential equation 



N 



V 2 V(r) = -s v e 2 J2 <K r i " *k) + s„e 2 n, j = l,2,...,N. 



(1.2) 




Here, the second term on the rhs comes from the particle-background interaction 
and n = iV/|A| is the mean number density of the mobile charges; the background- 
background interaction constant contributes to the partition function. 

The paper is organized as follows. In section 2, the method of the bulk Wigner- 
Kirkwood expansion is briefly reviewed in a format which can be relatively simply 
extended to the boundary case. The following sections are devoted to semi-infinite 
quantum systems constrained by a plane hard wall. In section 3, we propose a method 
for constructing the expansion of the Boltzmann density in configuration r-space. This 
expansion is subsequently used in section 4 for the generation of the expansion for 
statistical quantities, such as the partition function and the corresponding surface part 
of the free energy, one-body and two-body densities. Explicit calculations for Coulomb 
models with known classical statistical quantities are presented in section 5. A few 
concluding remarks are given in section 6. 

2. Wigner-Kirkwood expansion in infinite space 

Let the particle system be in the infinite space A = R u . In this so-called "bulk" regime, 
its equilibrium quantities in the nearly classical regime can be expanded in powers of h 2 
within the standard Wigner-Kirkwood expansion [U [2] which neglects fermion or boson 
exchange effects between quantum particles. In this section, we review shortly the 
derivation of this expansion by the Laplace transform method [H [16] which is suitable 
for our next purposes. 

2.1. Boltzmann density 

The Boltzmann density Bp in configuration space r can be formally written in the basis 
of plane waves as a i/A^-dimensional integral defined in an infinite domain R u : 



where p = (pi, p2, . . . , pjv) is the i/A^-dimensional momentum vector. Instead of 
considering e~@ H in the integration (12.11) . we take the Laplace transform of this operator 
with respect to the inverse temperature /?, 




(2.1) 




(2.2) 



and concentrate on the generation of an expansion in powers of % for 

„-(i/ft)p-r 1 _(i/ft)g-g 



(2.3) 
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Let us first rewrite H + z as 

H + z = D + 0, 
where D is a c-number 



D 



2m 



p 2 + V(r) + z 



and (9 an operator 
O -- 

and then expand 



— (-mvV 

2m V / 



2m 



P 



— - — e>— + — c— c— - 



F + 2 D D D D D D 
It can be verified that for any function /(r) the operator (9 acts as follows 

~ih _ - h 2 



o 



/(?)e 



(i/A)p-? 



_ e (i/ft)P'? 



m 



p v + 



2m 



One thus finds that 

-(i/fi)p-? 



,(i/R)p-? 



# + 2 

The term of order n in (12.91) 
1 



^ n=0 



ih 



m 



p-V 



n 2 -.2 
— v 

2m 



D 



D 



-p- V + — V 

m 2m 



D 



(2.4) 
(2.5) 

(2.6) 

(2.7) 

(2.8) 

(2.9) 

(2.10) 



is a polynomial in ft containing powers in the range from h n up to h 2n . This means that 
the truncation of the series (12.91) at some given order n provides all ft-terms up to order 
h n (or, equivalently, all A-terms up to order A n ). By carrying out all differentiations in 
(12,10]) . taking subsequently the inverse Laplace-transforms using the formula 

i r— iv'- 1 i 



Di (j - 1)! dzi- x z + p 2 /(2m) + V(r) 

1 



dpe 



-(3z_ 



-/3[p 2 /(2m)+F(r)] 



(2.11) 



(j = 1,2,...), and finally integrating on the momentum variables p, the Boltzmann 
density in configuration space (12.11) is obtained as the series 



re 



-m\ 



(2.12) 



n=0 



where B^\r) denotes the contribution implied by the nth term (I2.10p . In particular, 



Bf{r) 



1 

(V2^\)» N 
1 

{y/2n\) vN ' 



e-^A 2 
e-^|A 2 



(2.13a) 
(2.136) 



+ X 4 



24 1 



V 
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(2.13c) 



+ _w • V V) + -V (VV) + - V 

_|_ V y ■ v (w) - fj (w) vv + ^ (w) j | , 

etc. We conclude that the quantum Boltzmann density in configuration space is given, 
to order A 2 , by 



(r\e-P H \r) = 



1 



2ttA)^ 
1 



1 + A 2 



8-2 P z / ^ \ 



24 



+ 0(A 4 ) 



|e-' JV 






24 



A 2 -2 



24 



(2.14) 

When integrated over the coordinates of particle j, the corresponding term in 
J2f=i V 2 exp(— f3V) gives no volume contribution. Notice that odd powers of A do not 
appear in the Wigner-Kirkwood expansion, so that the bulk Boltzmann factor, when 
considered without the trivial prefactor (^/2tc\)~ uN , is invariant with respect to the 
transformation ih — > —ih. This is due to the fact that the integration of terms with odd 
powers of momentum components over the whole corresponding axis is 0. 



2.2. Statistical quantities 

According to the standard formalism of statistical quantum mechanics, the partition 
function of the iV-particle fluid (with ignored exchange effects) is given by the integration 
of the Boltzmann density over configuration space: 



l - ( dr (r|e-^|r 
V ! Ja x 



(2.15) 



For expressing macroscopic physical quantities, one defines the quantum average of a 
function /(?) as follows 



At the one-particle level, one introduces the particle density 



/ N 

M 1- ) = ( I>( r - r ^ 



\ j=l I qu 

At the two-particle level, the two-body density is given by 
< ) (r,r') = ( E 5{v-v 3 )5{v'-v k )\ . 

\ i,k=i / nu 



(2.16) 



(2.17) 



(2.18) 



It will be useful to consider also the truncated two-body density 

< )T (r, r') = n( 2 u )(r, r') - n qu (r> qu (r') (2.19) 
vanishing at asymptotically large distances |r — r'| — > oo, the pair distribution function 



£qu(r, r') 



n q ^(r, r / )/[n qu (r)n qu (r / )] and the pair correlation function h qu (r, r') = 
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<7 qu (r, r') — 1. The classical partition function and the classical average of a function 
/(f) are defined as follows 



Z = — 



dr 



N\ J a (V^rA)^ 



1 1 



df 



-/3V(r) 



-pv(f) 



/(?)• 



(2.20) 
(2.21) 



/ 27rA) l/iV 

According to our convention, the classical values of statistical quantities will be written 
without a subscript, like n(r), n( 2 )(r, r'), etc. 

Substituting the A-expansion of the Boltzmann density (I2.14p into formula (I2.15p . 
the quantum partition function takes the expansion form 

24* 



'qu 



1-A 2 ^(VV) + 0(A 4 ) 



(2.22) 



It is important to note that the general applicability of the Wigner-Kirkwood expansion 
scheme depends on the particular form of the interaction potential V; a divergence of 
the classical average (V V) indicates a non- analytic A 2 -expansion of the free energy. 

For the model system of the one-component plasma in v dimensions, the summation 
over the particle index j of the set of Poisson equations (11.21) leads to 



N 



V V(r) = -s u e 2 J2 S ( r j ~ r k) + s v Ne' 



n. 



(2.23) 



j,k=l 



The first term on the rhs of (I2.23|) does not give any contribution to (V V), because it 
is weighted by a classical Boltzmann factor which vanishes at zero interparticle distance. 
The replacement of V V in equation (12.22ft by the particles-background term s v Ne 2 n 
leads to 

r f3s v e 2 n 



1 - X 2 N- 



24 



+ 0(X 4 



(2.24) 



For the free energy F qu defined by the relation — /3F qu = lnZ qu , one gets from (12.241) the 
leading A 2 -correction term 



qu 



F + X N- 



s v e n 
24 



+ 0(A 4 



(2.25) 



which is, as one expects in the bulk regime, extensive in the number of particles N. 
Corrections to the free energy of higher order in A are standardly obtained from the 
partition function with higher-order A-terms by using the cumulant expansion (see, e.g., 
Ref. |!6j). 

The Wigner-Kirkwood formalism ignores the exchange effects. However, for the 
one-component plasma, they are exponentially small [17J. 

The A-expansion of the particle density (I2.17P for the one-component plasma reads 

A 2 



1 + ^V 2 + 0(A 4 ) 



rur 



(2.26) 
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In the bulk regime with the homogeneous classical density of particles n(r) = n, one 



obtains the trivial result n qu (r) 



n. 



The radial distribution function of the bulk one-component plasma has a A- 
expansion of the form [T8] 



X 2 

-V 2 

12 



0(X 4 



9{r) 



(2.27) 



In the short- distance r —>■ limit, g(r) is proportional to the Boltzmann factor of 
the bare Coulomb potential, which causes the breaking down of the Wigner-Kirkwood 
expansion scheme for too small values of r [18] . 



3. Half-space geometry: Boltzmann density 



Let the particles be constrained to the half-space A defined by Cartesian coordinates 
r = (x, r 1 ), where x > and r 1 - e R^ 1 denotes the set of [y — 1) unbounded 
spatial coordinates normal to x. The hard wall in half-space x < is considered to be 
impenetrable to particles. We group the components of the z/iV-dimensional vector in 
configuration space r = (r 1; r 2 , . . . , r N ) into two another vectors x = (xi, x 2 , ■ ■ ■ , x^) 

the corresponding momentum vectors generated from 

_L 



^); 



and i 

P = (Pi,P2, ■ ■ ■ ,Pn) are = (pf,p% , . . . ,p%) and p x = (p^, p^-, . . . , p^). The 
Hamiltonian is still of the form fll.ip . where the interaction potential ^(r) may include 
also additional effective one-body and two-body terms induced by the hard wall, e.g., 
the Coulomb energy of images in the case of charged particles near a dielectric wall. 

Due to the presence of the impenetrable hard wall, the wave function of the particle 
system has to vanish as soon as one of the particles lies at the wall. The one-body 
Boltzmann density in configuration space can be thus written in the basis of stationary 
plane waves with zero Dirichlet boundary condition at x = 0: 

dp* f dp 

/o 



(7Th) N J (27Th)(^ N 



re 



(3.1) 



where 



N 

n 



V2 



sin 



PjXj 



(3.2) 



To generate the ^-expansion of the matrix element (13.11) . we introduces the Laplace 
transform of the operator e _/3H and consider 

1 



H + z 



(3.3) 



One can proceed formally as in the previous bulk case, up to the expansion ( 12. 7(1 . 
Afterwards, one needs the following generalization of the operator equation (12.81) : 



O 



/(r> 



(i/A)p'-? 



,(i/?i)p'-? 



— p 

m 



v — v 2 

2m 



(3.4) 
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which is valid for any i/iV-dimensional momentum vector p' such that |p'| 
consequence, one finds from (12. 7p that 

-2 



iPl 



1 



H + z 



,(i/ft)p'-F 



-I oo 

s (i/»)r-?jL 



D 



E 

n=0 



ih -* h~ -2 
-p' • v + —V 

m 2m 



As a 



(3.5) 



To make use of this operator formula, we consider an exponential representation of the 
function i/jp(r) given by equation (13.21) : 



(x/2i 



N 



N 

n e 

i=i CTj =±i 



^■e 



(V2i 



A' 



En 



,(i/a)p(a)-r 



(3.6) 



{<?} 



Here, each of the cr-components of the "state" vector a = (a\, o%, . . . , a^) can take 
one of the values ±1 (the sign determines the direction of the plane wave propagation) 
and the summation goes over all possible 2^ a-configurations. For a given state vector 
a, we have redefined the particle momenta Pj(cj) = (ajPj,pf) (j = 1,2, ...,N) and 
grouped them into the z/iV-dimensional vector p(cr) = (pi(o"i), p 2 (o"2), . . . , Pn(o~n)), 
which couples via the scalar product to the r-vector in the exponential of equation 
(13. 6p . Since the momentum vector p(<?) satisfies the requirement |p(er)| = |p| for an 
arbitrary state vector a, one can apply the key relation (13.51) to every summand in the 
exponential representation ( 13. 6p of ^(r) to obtain 



pv 'H + z 



^(?) = ^(F) 



N 



X 



D 



(V2i 

oo ( 

E 

n=0 I 



N 



En 



s (i/ft)p(aO-? 



{<?} 

ih^,_>. -l h 2 -2 

-p(<x)- V + — V 

m 2m 



Let us consider the term of order n in the expansion (13.71) 



N 

n 







sin ^ 






first 


carries 




P( 



1 w 



{*} j=l 




ih 



m 



-p{a) 



V + — v 2 

2m 



One first carries out all differentiations in (13. 8p . Since the scalar product 







N 



E^ + Ep^-v, 



(3.7) 



(3.J 



(3.9) 



fc=i 



contains components of the state vector a, each of the components can take either odd or 
even power multiplied by the corresponding exponential; the summation over o~j = ±1 
then gives rise to the original stationary wave oc sm{p x -Xj/h) for odd powers of a 3 - and 
to the shifted wave oc cos(pjXj/h) for even powers of o~j. As the next step, one performs 
the inverse Laplace transform of (j = 1,2,...) using the formula (12.111) . In the 

consequent relatively simple integrations over momenta we have to keep in mind that 
while the integration of odd powers of the p x components give zero contributions, this 
is no longer true for odd powers of p x since the integration goes only over the positive 
values of p x . The final result for the Boltzmann density in configuration space (13. ip has 
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the form analogous to that of the bulk case ( 12. 12}) . (r|e 1311 |r) 



(n), 



r) where 



.Bg^(r) denotes the contribution implied by the nth term (13.81) . 



P 



We have performed the outlined procedure for the first terms B^\ Bp and B^\ 
We do not go into details of the calculations, but only write down the final result and 
then discuss the origin of the terms obtained and their structure. 



The result for B 



(0). 

p 1 



r) reads 
dp* 



dp 



_L 



o (7rh) N J {2nh)(»-V N j 



N 

n 



2 sin' 



PjXj 



-/3[p 2 /(2m)+V(r)} 



1 



-pv 



N 

n(i 



(3.10) 



-2* 2 /A 2 \ 

^ttA)^' >' 

In comparison with the corresponding bulk term (12.13aj) . each particle gets an additional 
"boundary" factor 1 — exp(— 2x 2 /X 2 ) which goes from at the boundary x — to 1 
in the bulk interior x — > oo on the length scale ~ A. The product of boundary factors 
ensures that the quantum Boltzmann density vanishes as soon as one of the particles lies 
on the boundary. The dependence of the boundary factor on the de Broglie wavelength 
A is non-analytic; this fact prevents one from a simple classification of contributions 
to the Boltzmann density according to integer powers of A like it was in the bulk case. 
However, as will be shown explicitly in the next section on the evaluation of the partition 
function, when the exponential part exp(— 2x 2 /A 2 ) of the boundary factor (multiplied 
eventually by another A-independent function of x) is integrated over the x-coordinate, 
the analyticity in the parameter A is restored. At this stage we only notice that the 
integrated product of boundary factors should be expanded in the following way 

N N i N 



3=1 



nCl-e-^'^l-Ee-^ + i £ 6-^6-^ + - ••■(3.11) 



i,fe=i 



where each exponential term exp(— 2x 2 /A 2 ), when integrated over x, produces one A- 
factor as the result of the evoqued substitution of variables x = Xx'. 



The result for B 



P 1 



B 



r) is found in the form 
1 



N N 



(V^fA)^ 

N 

+n(i- 



-pv. 



-2x 2 /A i 



e n (i 



-2x 2 /\ 2 







k=l 3=1 



dV 

dx k 



vv 



(3.12) 



Here, the first term corresponds to the source operator D _1 p(cr) ■ VD _1 whose bulk 
counterpart does not contribute to the Boltzmann density (12.1361) . The dependence 
on A appears, besides the product of the particle boundary factors 1 — exp(— 2x 2 /X 2 ), 
also via the combination x exp(— 2x 2 /X 2 ). This function has a maximum of order A 
and the integration of this function over x gives a contribution of order A 2 . These 
properties make the function xexp(— 2x 2 /X 2 ) of an order equal to or "weaker" than A, 
and therefore this function is a legitimate expansion parameter. The second term on the 
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^2 . 



rhs of equation f l 3 . X 2 j) . having the origin in the source operator D _1 V D _1 , corresponds 
to the simple multiplication of its bulk counterpart (12.1361) by the product of particle 
boundary factors 1 — exp(— 2x 2 /\ 2 ). 



The final formula for BfP(r) consists in four terms 



Bf{v) = B i 



'3 

(2,1), 



B?*m + B™m + B™m 



(3.13) 



where 



Bp ; (r) 



-J e-^( TT (l - e- 2 ^) A 2 [^vV - (VV) 2 



e •? 



N N 

+e n (i 

k=l i=i 
CM*) 

N N 

e n (i - e "" ?/A2 

fc,!=l 3=1 
(MO CMM) 



2 * 2/A ^ x 2 e- 2 "^ A2 



2/? d 2 T/ _ /J 2 / ay ' 

3 cte? 2 Ictefc 



x fc e- 2 ^-/ A2 x ie - 2 ^/ A2 



2/5 «9 2 K 



(3.14a) 
f3 2 dV dV 



3 dxkdxi 2 ctefc <9x; 



corresponds to the source operator D 1 p(?) ■ V[D 1 p(o ? ) ■ VD 1 ], 



<)(r) 



■i AT N 

= e-^V TT fl-e- 2 ^ 2/A2 )x fc e- 2 ^/ A2 A 2 



fc = l 3 = 1 

(M*0 



/3 a 

6 dxi 



+^r\ V~\' 



w + CA(w) 2 -^(v^ ! " 1 



Sxi. 12 ctei 



^2 



(3.146) 



has its origin in the source operator D 1 p(&) • V[D X V -D 



-ii 



< 3 »(?) 



AT AT 



^A) 



+^ ( vV) — + f -— (vvY - f - fw 



uN 



e n (i- 

k=l 3=1 
(Mfe) 

9^ /3 2 5 



-2x 2 /A : 

e j 



-2a;2/A2 , 2 



P__d_ 

6 (9xi 



/? 3 



<9x 



8 dxi 



10 



2 

dx k 



^2 



(3.14c) 



comes from the source operator D 1 V [D 1 p(c?) • V-D 1 ], and 



< 4, (F) 



A 



-/3V 



2ttA) 



i/JV 



I] 1 - e"^ 2/A2 A 



£ 2 



16 

/3 3 



+ > w -v v v) + —v w +— vv 



18 



24 



£ 2 



32 



_^_w ■ v (w) z - 1 (w) z vV + £ (vy 



2 ^2 



30 



24 



72 



(3.14d) 



results from the source operator D X V [D *V D 1 ]. It is easy to check that the term 
Bp(r) as a whole is of the leading order A 2 , like it was in the bulk regime (see equation 

pa ). 

We would like to note that in the bulk interior far away from the boundary, i.e., when 
all x-coordinates of particles {xj — > oo}^ =1 , the results for the boundary Boltzmann 
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density must reduce to the bulk ones fl2.13ap - fl2.13cp . It is easy to verify that the 
results obtained pass this test. 

To summarize, the A-expansion of the Boltzmann density in the presence of 
a boundary is more complex than the analytic one in the bulk. The expansion 
involves not only powers of A, but also non-analytic position-dependent terms of type 
1— exp(— 2s 2 /A 2 ), a;exp(— 2x 2 /A 2 ), etc. The first three Boltzmann terms B^\f), B^\r) 
and B^\r) exhibit properties similar to their bulk counterparts ( I2.13ap -( 12.136P : The 
maximum of B^\r) is of order A™. We anticipate that this formal structure of the 
Boltzmann terms is also maintained on higher levels. 



4. Half-space geometry: statistical quantities 

4-1. Partition function and free energy 



Since the Boltzmann density is represented by the series 

E~=o^ n) (?), the partition 



function (I2.15P is expressible as 

^u = E4?> Z^ = — dvBf\v). (4.1) 

n=0 ly ■ JA 

In the next paragraphs we evaluate consecutively the A-expansion of the first three terms 
Zj®, ZjV and Z^ , up to the order A 2 . According to the analysis in the previous section, 
higher-order Z$ terms are expected to contribute only to A 3 and higher powers of A. 

Inserting the Boltzmann term B^ f l3.1Up into the definition (14. ip of Z$® and 
performing the expansion (13. lip for the product of boundary factors, we obtain 

Z2L = 1- dre- 2 * 2 / A2 n(r) + -i / dr x / dr 2 e-^V^V^ri, r 2 ) + • • ■ . (4.2) 
Z J a 2! Ja Ja 

The system is translationally invariant in the v 1 - space, i.e. n(r) = n(x), n.^(ri, r 2 ) = 
n( 2 \x\,X2] \r± — r 2 |), etc. This property enables us to rewrite equation (14. 2 j) in a more 
convenient form 

roo , , , 1 r poo 

= 1 - I9AI / dxe~ 2x /A nix) + -=-|c>A| 2 / dx< 
Z Jo 2! [Jo 

I P P OO P OO 

+ ^\dA\ dr ± dxi / dx 2 e- 2a;2/A2 e- 2 ^/ A2 n (2)T (x 1 ,x 2 ;|r ± |) + ---, (4.3) 
2! J Jo Jo 

where \dA\ denotes the surface of the half-space domain A at x = 0, and the third 

term has been added to the rhs and subsequently subtracted from the double integral 

(inducing in this way the truncation of the two-body density) in order to arrange the 

series for the application of the cumulant method. Assuming that the classical density 

profile n(x) is analytic at the boundary x = 0, i.e. n(x) = n(0) +n' (0)x+n" (0)x 2 /2! + - • • 

(this condition is not always fulfilled, see subsection 5.1.), the integral 

1 „ , „ poo , 

dxe~ 2x /A n{x) = A / dx'e" 2 * n(Ax') 



e ~ 2x2/x2 n(x) 



A^/|M0) + A 2 ^(0) + O(A 3 ). (4.4) 
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Performing an analogous procedure in the double integral on the rhs of equation (14.3)) . 
we finally arrive at 



qu 



l-|0A|Aiyfn(O) + I 



\dA\X^n(0) 



|<9A|A 2 -n'(0) 



4 



+ \dA\X 2 ^Jdr ± n^ T (0,0;\r ± \)+O(X 3 



(4.5) 



We recall that the boundary values of the classical statistical quantities n(0), n'(0) and 
n^ T (0, 0; \^ ± \) are, in general, nonzero. 

The boundary factors 1 — exp(— 2x 2 /X 2 ) do not play any role in the Boltzmann 
term Bf (l3~T2l when one is interested in the A-expansion of the corresponding Zf^ up 
to the A 2 order, so we can ignore them. After simple algebra, Z$ is written as 



qu 



r dre -^ 2 /* 2 xn >( x ) + a 2 

Ja 



•f(vV)4((wf 



+ 0(A 3 ), 



(4.6) 



where the classical average is taken over the half-space A. The integral on the rhs of 
( 14. 6 p can be treated in the way outlined in the previous paragraph, with the result 

,1 



J dre- 2x2 / x2 xn'(x) = |<9A|A 2 -n'(0) + 0(A 3 ), 
the second term can be simplified by using the evident equality 
f3 2 ({VV) 2 ) =/3(vV\+ / drV 2 n(r) 



(4.7) 



= (3{V - |9A|n'(0), 
provided that n'(0) is finite, see subsection 5.1. 



(4- 



From the four contributions to the Boltzmann term Bp (r) in (13.131) . only the first 



one implies the A term: 



■^qu 



A 



^2 

V V 



2 



yvy 



0(A 3 ). 



(4.9) 



Also here we can apply the relation (14.81) to eliminate the term (3 2 ((VV) 2 ) in favour of 
the Laplacian term. 

The results obtained in this subsection can be summarized by the expansion formula 
for the free energy, expanding in powers of A the logarithm in f3F c 



qu 



\nZ, 



qu- 



24 



PF qu = (3F + A 2 ^ < V'V ) + \dk\ \ X 



2 V 2 



n(0) 



+X 2 ±n>(0) - A 2 ^ / dr^ 2 > T (0,0; |r^|)} + 0(A 3 



(4.10) 



The free energy of a quantum particle system constrained to a domain A of volume |A| 
and surface |9A| possesses the following general form 

F qn = f b \A\+ lqu \dAl (4.11) 



where f^ u is the bulk free energy per unit volume and 7 qu is the surface tension. In 
the formula (14. lOj) . the quantum contribution to the bulk part of the free energy comes 



Wigner-Kirkwood expansion for semi-infinite quantum fluids 14 

-> 2 

only from the potential-dependent term oc (V V). This term can, in principle, give also 
some surface contribution. But for our model system of the semi-infinite one-component 
plasma constrained by the plain hard wall, the potential V fulfils the Poisson equation 
(12.23P and therefore the term oc ( V V) produces only the bulk quantum contribution to 
the classical free energy, as shown in equation (12.251) . By comparing (I4.10p with (14. lip , 
the surface tension is found to be 

pj q u =pi+ 4/!^°) + x2 k n ' {0) - x2j k I dr± ™ (2)T(o ' 0; |r±|) + o(A3) - (4 - i2) 

It is interesting that up to the A 2 order the interaction potential V enters into the 
formula for the surface tension only implicitly via the classical one-body and two-body 
averages. We note that, although for the one-component plasma the boundary truncated 
two-body density n^ 2 ^ T (0, 0; | r — 1 ) in (14. 121) exhibits a long-range decay of type l/|r J -| iy 
at asymptotically large | r- 1 - 1 — > oo [T51 [2D], this makes no divergence problem at large 
distances in the integration over the {y — l)-dimensional r 1 - space. 

Without going into details, the result for the surface tension (I4.12p can be 
straightforwardly generalized to many-component quantum fluids composed of different 
species a = 1,2, with masses {m a }^ =1 and the corresponding de Broglie 

wavelengths A a = K{(3 /m a ) l l 2 : 

Pl qu = Pi + A «^(°) + ^ E A X(0) 

" Y^/ dr± E A «ViT(aO; \v L \) + 0(A Q A /3 A 7 ). (4.13) 

a, (3 

We note that in the special case of a multi-component Coulomb fluid containing at 
least two species with charges of opposite signs, a short-distance regularization of the 
Coulomb potential is needed in the evaluation of the classical averages in (I4.13P in order 
to prevent the thermodynamic collapse. 



4-2. Particle density profile 

As concerns the quantum corrections to the classical density profile and two-body 
density, since they involve a quantum Boltzmann factor which is not integrated on all 
variables, their expansion will contain, in addition to powers of A, non-analytic terms 
exp(— 2x 2 /X 2 ). Here, we shall restrict ourselves to the leading correction of order A, 
considering a factor xexp(— 2x 2 /A 2 ) as being of order A since its maximum is of order 
A, as pointed out in section 3, using the expansion technique presented in the previous 
subsection 4.1. 

The one-body density profile (I2.17P is obtained in the form 



n qu {x) 



1 - e- 2 * 2 ^ 



nix 



-4/l/ dr±n(2)T(x ' 0;|r±|) 



the second term comes from Bf (?) (M) . 



+ xe~ 2x /x V(x) + 0(A 2 ), (4.14) 

.(o 

P 



where the first (product) term has its origin in the Boltzmann term Bg(r) (I3.10p and 



Wigner-Kirkwood expansion for semi-infinite quantum fluids 15 

In the classical description of the one-component Coulomb plasma, the system is 
neutral whole, 

dx [n(x) -n] = 0, (4-15) 



o 



and a particle of charge e is surrounded by a screening cloud the average charge of which 
is exactly — e, 



n[x) 



-J dr ± J°° dx'nW T (x,x'; (r^). (4.16) 

The neutrality sum rule (I4.15P is valid also in the quantum regime. For our purposes it 
can be reexpressed in the form 

dx [n qu (x) - n(x)] = 0. (4.17) 



It is easy to check that for the density profile (14.141) this relation is fulfilled in the linear 
A order. 

The formula (14.141) can be straightforwardly extended to a quantum fluid with many 
components a = 1, 2, . . . , L: 



, v 1 (W 
n a (x) - 2 



+ xe- 2x2/x «n' a (x) + 0(A 2 ). (4.18) 



4-3. Two-body density 



The leading quantum correction for the truncated two-body density fl 2 . 1 9 [) is found to 
be 

ng} T (n,r 2 ) = fl -e- 2 ^ A2 ) (l-e-^ x2 )\n^ T {v v r 2l 



'(n,r 2 )= (l-e- 2 ^ A2 )(l-e- 2 ^ A2 )|, 

-4/l/ dr±n(3)T[ri ' r2 ' (o ' r±)] 



1 _ e -2*|/A 2 A a . ie -2^/A^ n (2)T (ri; , (419ffl) 
' OX\ 

1 - e" 2 ^ 2 ) x 2 e"^/ A2 A n (2)T (ri; } + 0(A 2 } 
' ox 2 



where 



n (3)T (n, r 2 , r 3 ) = n (3) (n, r 2 , r 3 ) - n (2)T (n, r 2 ) n(r 3 ) - n (2)T (r 1; r 3 ) n(r 2 ) 



- n (2)T (r 2 , r 3 ) n(r x ) - n(ri)rc(r 2 )n(r 3 ). (4.196) 



5. Coulomb models with known classical statistical quantities 

5.1. High-temperature Debye-Hiickel limit 

For the semi-infinite one-component plasma, at least in two and three dimensions, the 
classical truncated two-body density and density profile are known in this limit (see, 
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e.g., equations (3.12) and (3.27) in [12]). The classical surface tension is also known (see 
equations (4.6) and (4.7) in [21]). The integrals which appear in (I4.12p and (I4.14p are 
easily computed: 

J drV 2)T (x, 0; (r^l) = -nKe~ KX , (5.1) 

where k = [2(u — l)7r/3e 2 n] 1/ ' 2 is the inverse Debye length. It should be noted that the 
density profile 

n(x) = n + 




r e— - 2(1 + t 2 ) 1 / 2 exp[-2(l + t 2 ) 1 / 2 ^] 2 

Jo [(l + t*)W+t)*(3 + 4t*) dt ' ^ 

contrarily to the assumption in subsection 4.1., can be expanded in powers of x only up 
to order x 3 ~ u : the nth derivative of n(x) is infinite for n > 4 — v. 

In two dimensions, n'(0) is finite and the quantum surface tension up to order A 2 
( I4.12p becomes 

„_f ln3 _ 1 + wiW 

V 9 / 8tt 

+ A 2 (^ + ^)+o(A 2 ). (5.3) 

It is likely that ^"(O), which is infinite, would appear if we attempted to compute the 
term of order A 3 . This is an indication that the Wigner-Kirkwood expansion of the 
surface tension is not in integer powers of A in the present case; some singular term 
(logarithmic for instance) appears beyond the order A 2 . 

In three dimensions, n'(0) is infinite. The integral (14.41) is nevertheless finite, 
although it cannot be expanded up to order A 2 (instead, a term of order A 2 In A appears). 
Also, using (14. 8p in the contributions (14.61) and (14. 9p to the partition function shows 
that these contributions, formally of order A 2 , are infinite. This is an indication that 
some singular term, which we are unable to estimate, appears in the Wigner-Kirkwood 
expansion of the surface tension beyond the order A. At order A, (I4.12p gives 

21n2 - 1 o . 1 [W 



Plau = ft 2 + A- 

f Iq 32tt 2 V 2 



fl - 31n3 + 7r\/3 k 3 
n — 



24vr 



+ o(A). (5.4) 



For the two-component plasma, the classical density profiles (see equation (4.9) in 
[j~9]). which contain the singular Bessel function K2(2kx), also cannot be expanded in 
powers of x to all orders. 

It is also interesting to compare the density profiles (14.181) for a many-component 
semi-infinite fluid with the results of Aqua and Cornu [T2] , which have been obtained by 
a very different method, for a three-dimensional multi-component semi-infinite quantum 
Coulomb fluid. Let the charge of a particle of species a be e Q . The inverse Debye length 
now is k = (4:iiPJ2a n a e a) 1 ^ 2 - An easy generalization of (15. ip shows that 

dr^g T (x, 0; Ir-H) = _ 4 ^"Wy ^ (5 5) 

K 



I 
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Thus, ( 14.181) becomes 



ng?(*)= (l-e-^) 



\/27r 3 / 2 /5n a e a 
w a (a?) H 2^ A 



7 



+ xe" 2:!; /A V Q (x) + 0(A 2 ). (5.6) 

The first line of ( 15. 6ft is identical to the result (1.14) in Aqua and Cornu [12]. However, 
they do not have the second line, since this second line is of higher order in nf3e 2 and 
k\ than in the regime they have considered. More generally, terms of higher order in A 
are also of higher order in the classical dimensionless coupling parameter nf3e 2 and, in 
this weak coupling regime considered by Aqua and Cornu, there are no terms o(A) in the 
density profiles. 

5.2. Exactly solvable two-dimensional one- component plasma 

The classical two-dimensional one component plasma is an exactly solvable model when 
(3e 2 = 2. In particular, for the semi-infinite geometry, the density profile and the two- 
body densities are exactly known [19]. Now, the density profile 



2 exp - (t - xV2 2 1 
n(x) = n^= / FL / ' J dt, 5.7 

v ' Jo l + $(t) V ; 

where $ is the error function, can be expanded in integer powers of x: all the derivatives 
of n(x) at x = are finite (this might be a property of even integer values of (3e 2 = 2). 
In particular, 

n(0) =n In 2, (5.8) 
n'(0) = - v^(2n) 3 / 2 jH In (^f^j dt, (5.9) 

where the integral has the numerical value -0.3377. The integral which appears in ( 14.121) 
can be computed: 



/d r V^(0,0;| r -|)^-4^ n 3 /2 ^_^_ dti (5 , 0) 

where the integral on t is related to the previous one: by integrations per partes, one 
finds 

e" 2 * 2 , ^ ^ f» /l + $(t)\ , 

dt = ^— + - In v ; dt = 0.3558. (5.11) 



lo [l + $(t)] 2 2 2 Jo V 2 

The classical surface tension is also known |21j : 

^=_j]5:r ta fi±*wu. (5.12, 



2tt Jo V 2 

Using (I5TT21 . (15131) . flES]) , and (15101) in (147L21) gives the surface tension. 
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6. Conclusion 

For nearly classical quantum fluids (nearly classical means that the quantum effects are 
weak), the Wigner-Kirkwood formalism can be generalized to semi- infinite fluids, i.e. 
to fluids confined to a half-space x > by a hard plane wall at x = 0: the boundary 
condition for the wave functions is that they have to vanish at x = 0. Under the 
condition that the classical density profile has all its derivatives with respect to x finite 
at x — 0, the surface tension can be expanded in powers, including odd powers, of 
the thermal de Broglie wavelength A; this is to be compared with the case of infinite 
homogeneous fluids with sufficiently smooth interactions, which have a free energy 
the expansion of which contains only even powers of A. The one-body and many- 
body densities are more complicated. Their expansions contain not only powers of A, 
but also non-analytic position-dependent terms, localized near the wall, of the types 
exp(— 2x 2 /X 2 ), xexp(— 2x 2 /A 2 ), etc. 

The above assumption of the classical density profile having all its derivatives with 
respect to x finite at x = is not always fulfilled; in particular, it fails for Coulomb 
fluids in the high-temperature Debye-Hiickel regime. For such cases, the surface tension 
can be expanded in integer powers of A only up to certain order, then singular terms 
appear, for instance of the type A n In A. 

As a task for the future, one should consider the exchange effects for the half-space 
system geometry. We anticipate that, for the considered one-component plasma, the 
exchange effects are exponentially small and therefore negligible in the nearly classical 
regime, like in the bulk case [PT] . 

The extension of the present formalism to more general boundary conditions is 
needed. For example, in the case of the Neumann boundary condition at x = one has 
to write down the Boltzmann density in the basis of waves \/2cos(p x x/h) exp(ip ± -r ± //i) 
which derivation with respect to x vanishes at x = 0. Other geometries of the confining 
domain (e.g., a strip of finite width) are also of interest. 

As concerns the model of the one-component plasma, we plan to extend the study to 
the hard walls made of conducting and dielectric materials; such walls imply additional 
image forces acting on charged particles. In connection with Coulomb fluids, one should 
also check explicitly all available quantum sum rules for the half-space geometry. 
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